Spatially Resolved Millimeter Interferometry of 
SMM J02399-0136: a Very Massive Galaxy at z = 2.8 
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ABSTRACT 



We report high-resolution millimeter mapping with the IRAM Plateau de 
Bure interferometer of rest-frame 335 /xm continuum and CO (3-2) line emission 
from the z = 2.8 submillimeter galaxy SMM J02399-0136. The continuum emis- 
^ ■ sion comes from a ~ 3" diameter structure whose elongation is approximately 

east-west and whose centroid is coincident within the astrometric errors with the 
>- ■ brightest X-ray and rest-UV peak (LI). The line data show that this structure is 

^ ■ most likely a rapidly rotating disk. Its rotation velocity of > 420kms~^ implies 

' a total dynamical mass of > 3 x 10^^ sin~^ i Hqj Mq within an intrinsic radius of 

^ 8 h^^ kpc, most of which is plausibly in the form of stars and gas. SMM J02399- 

g ■ 0136 is thus a very massive system, whose formation at 2; ~ 3 is not easy to 

■ understand in current CDM hierarchical merger cosmogonies. 

^ \ Subject headings: galaxies: formation, kinematics and dynamics, active, star- 

V-j \ burst; cosmology: observations 

c/3 ' 



1. Introduction 

The strength of the extragalactic far-IR/submillimeter background indicates that about 
half of the cosmic energy density comes from distant, dusty starbursts and AGN (Fixsen 
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et al. 1998; Hauser et al. 1998; Lagache et al. 1999; Pei, Fall, & Hauser 1999). Surveys 

with ISOCAM at 15 /xm, SCUBA at 850 /xm, and MAMBO at 1200 fxm suggest that this 
background is dominated by luminous and ultraluminous infrared galaxies ([U]LlRGs: Lin ~ 
IQii-s _ 10^3 Lq) at ^ > 1 (e.g., Small, Ivison, & Blain 1997; Hughes et al. 1998; Barger, 
Cowie, & Sanders 1999; Aussel et al. 1999, 2000; Carilh et al. 2001; Elbaz et al. 2002; 
Genzel & Cesarsky 2000). Faint submillimeter sources frequently have only weak (if any) 
counterparts in the rest-frame UV and optical, however (Small et al. 2000, 2002; Dannerbauer 
et al. 2002). Redshifts and spectroscopic parameters have thus far been confirmed with CO 
interferometry for only three of the > 10^ detected systems (Blain et al. 2002, and references 
therein) . 

One of these three sources is SMM J02399-0136 (hereafter J02399), the brightest of the 
15 significant background sources discovered in the SCUBA cluster lens survey (SCLS: Small 
et al. 1997, 2002). This galaxy's 850 //m flux density of 23.0 mJy includes a magnification 
factor of 2.45 due to the z = 0.37 foreground cluster Abell370 (Ivison et al. 1998; R. Ivison, 
private communication). For a dust temperature of 45 K and emissivity index f3 = 1.5 
(typical for local luminous IRAS galaxies: Dunne et al. 2000), its intrinsic IR luminosity 
(from 8 — 1000 /um in the rest frame) is Ljr = 1.2 x 10^^ Lq.^ J02399 is rich in gas as 
well as dust: Frayer et al. (1998) detect strong CO(3~2) emission at zco = 2.808, slightly 
redward of its UV hne emission at zyy — 2.803 (Ivison et al. 1998). The optical through 
X-ray properties of J02399 resemble those of local type 2 AGN and broad absorption hne 
QSOs (Ivison et al. 1998; Vernet & Cimatti 2001). Moderately broad (2000 km s"^) Lya 
emission is extended over about 13", with a maximum on the brightest UV continuum peak 
LI, and a second maximum on the UV continuum peak L2 lying ~ 3" farther east (Ivison 
et al. 1998; Vernet & Cimatti 2001). Bautz et al. (2000) find that strong X-ray emission 
from a compact nucleus is obscured by a column of ~ 10^^ cm^^ and estimate an unabsorbed 
X-ray luminosity of 2.7 ± 1.0 x 10^^ h^'\ erg s~^. The ratios of X-ray, far-IR, and radio fiuxes 
suggest that 15-70% of the bolometric (infrared) luminosity may come from the AGN (Bautz 
et al. 2000; Frayer et al. 1998), as is true for composite starburst/AGN members of the local 
ULIRG population hke Mkn231 (Sanders & Mirabel 1996; Genzel et al. 1998). 



^This paper assumes a flat Oa = 0.7 cosmology with //q = 70 /lo.ykms ^ Mpc ^, for which zq,o = 2.808 
corresponds to a luminosity distance £>l = 23.5 K^Jj Gpc and an angular diameter distance Z)a = 1-6 /i^y Gpc. 
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2. Observations 

We observed J02399 with the B, C, and D configurations of the IRAM Plateau de Bure 
interferometer (PdBI: Guilloteau et aL 1992) in November - December 1998 and in Jan- 
uary - September 2002. In 1998, the array consisted of five 15 m telescopes, each equipped 
with both [single-sideband] 3 mm and [double-sideband] 1 mm SIS receivers which we used 
simultaneously. System temperatures (referred to above the atmosphere) were ~ 120 K and 
~ 400 K at 3 mm and 1mm, respectively. We tuned the 3 mm receivers to 90.81GHz, the 
frequency of the redshifted CO(3-2) line for zco — 2.808, and the 1 mm receivers to a 
line-free region at 235.00 GHz corresponding to a rest wavelength of 335 /im. Four correla- 
tor modules were deployed to sample the full 536 MHz of contiguous 3 mm bandwidth at 
2.5 MHz channel"^ resolution, while the remaining two modules provided 295 MHz of 1 mm 
continuum bandwidth. By 2002, the array had expanded to six telescopes, eight correlator 
modules, and a wider IF bandwidth. We therefore retuned the 1 mm receivers to 211.83 GHz, 
the frequency of the redshifted GO(7-6) fine, and deployed four correlator modules (giving 
559 MHz of continuous bandwidth) at both 3 mm and 1 mm. The pointing center was also 
changed by 0.91" for the 2002 observations, to agree with the GO(3-2) centroid {a (J2000) 
= 02'*39"^51.89^ and S (J2000) = -01°35'59.9") measured by Frayer et al. (1998). 

We calibrated the 1998 and 2002 data separately using the GLIG routines in the IRAM 
GILDAS package (Guilloteau & Lucas 2000). Passband calibration used one or more bright 
quasars. Phase and amplitude variations within each track were calibrated out by interleav- 
ing reference observations of one or two quasars near the source every 30 minutes. The overall 
fiux scale for each epoch was set by comparisons of model sources CRL618 and MWC349 
with bright quasars whose flux densities are regularly monitored at both the PdBI and the 
IRAM 30 m telescope. Once the source data were calibrated, we precessed the 1998 visibil- 
ities to the 2002 phase center. After confirming our nondetection of CO (7-6) emission, we 
rescaled the 2002 1mm continuum visibilities from 211.83 to 235.00 GHz using a power law 
index of 3 derived from (sub)miUimeter photometry (Ivison et al. 1998; Small et al. 2002). 
We then combined the 1998 and 2002 uv data, yielding a total of 40 distinct basehnes (24- 
331m) from the three array configurations, and smoothed the CO(3-2) data to 15kms~^ 
resolution. After fiagging bad and high phase noise data, we were left with the [on-source, 
six-antenna array] equivalents of 39.6 hours of 3 mm data and 21.6 hours of 1 mm data. We 
deconvolved the data cubes with natural weighting using the CLEAN algorithm as imple- 
mented in the NRAO AIPS package (van Moorsel, Kemball, & Greisen 1996) and GILDAS. 
Although no 3 mm (rest-frame 840 /^.m) continuum emission was seen, both the CO (3-2) fine 
and the [rest-frame] 335 /im continuum emission are well detected and resolved in our maps. 
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3. Results 

3.1. Continuum emission 

Figure 1 shows the 335 //m (Aobs — 1.27 mm) continuum emission from J02399 at a 
resolution of 1.8" x 1.4", superposed on a sub-arcsecond resolution it!-band image obtained 
at Keck (L. Cowie, private communication). The observed total flux density of 7.0±1.2mJy 
(±1.4 mJy for a 20% uncertainty in the flnx scale) agrees well with the value (8.5 mJy) 
expected from JCMT photometry at 850 and 1350 fira, suggesting that the interferometer 
has resolved out little if any emission. The central continuum structure is clearly resolved in 
the east-west direction, parallel to the dominant lensing shear predicted by the SCLS lensing 
model (R. Ivison, private communication). A fit to the uv visibiUties yields FWHM diameters 
of 3.2" (east- west) and < 0.8" (north-south); correcting the former for magnification by a 
factor oi M. — 2.45 imphes a source-plane hnear size of ~ 10/i^ykpc. An additional tail of 
emission extending > 4" WSW from the central peak may be associated with extended Ly a 
emission detected by Vernet & Cimatti (2001), and with another fainter i?-band emission 
peak (Figure 1). The ~ lO/igykpc size of the central dust source in J02399 is about three 
times larger than the sizes of local ULlRGs (Downes & Solomon 1998; Bryant & Scoville 
1999). This fact may not be surprising, however, given that J02399 is about ten times 
more luminous than the typical local ULIRG. Since the star-forming interstellar medium in 
ULIRGs is thought to consist of a high volume filhng factor distribution of dense gas and 
dust (Downes & Solomon 1998; Sakamoto et al. 1999), a ten times larger star formation rate 
requires a ten times larger volume of star-forming regions. 

Table 1 summarizes the absolute positions and astrometric errors of the emission peak 
of J02399 in the different wavebands. The coincidence of the far-IR centroid with the UV 
and X-ray continuum peaks within the combined astrometric uncertainties of ~ 1" suggests 
that LI is indeed the primary origin of J02399's dust emission. The (bluer) UV continuum 
and LycK emission peak L2, in contrast, has no associated dust emission; this supports the 
idea that L2 is primarily scattered AGN light from LI (Vernet & Cimatti 2001), rather than 
a second galaxy in the early stages of a merger with LI (Ivison et al. 1998). In fact, there 
may be an anticorrelation between the fainter rest UV emission and the 335 fim continuum 
emission, suggesting that extinction strongly affects the short-wavelength morphology of 
the source. UV emission may only escape from "extinction holes" in the extremely dusty 
environment of J02399, a scenario also suggested by Ivison et al. (2001) in their study of the 
dusty high-redshift SCUBA source SMM J14011+0252 (hereafter J14011). 



- 5 - 



3.2. CO line emission 

Figure 2 shows the velocity-integrated CO(3-2) hne emission (Aobs = 3.30 mm) at a 
resolution of 5.2" x 2.4" in greyscale, with the continuum contours of Figure 1 overlaid. 
Taking into account the different spatial resolutions, the distributions and centroids of the 
velocity-integrated CO (3-2) and 335 /xm continuum emission are coincident, and the peak 
CO emission (like the peak continuum) originates in the eastern (blue-shifted) emission peak. 
The UV/X-ray source LI is located 0.7" ± 0.6" from the rotation center. Figure 3 shows the 
CO (3-2) line profile, smoothed to 60kms^^ and integrated over the central ~ 5". Figure 
4 displays channel maps of the CO(3-2) emission (spaced by 200kms~^) at a resolution 
of 5.2" X 2.4". Our new data agree with and significantly improve on the earlier 7" x 5" 
interferometry of Frayer et al. (1998). The molecular line emission of J02399 is very broad 
and has a symmetric "double-horn" profile which is characteristic of rotation. Line emission 
extends over more than 1100 km s~^ (FWHM), and the global profile shows emission peaks 
at — 520kms~^ and -1-230 km s~^ (relative to the assumed zco — 2.808). An improved CO 
systemic redshift for the source would thus be 2.8076 ± 0.0002. Channel maps show that 
blue and red-shifted emission maxima originate 1.5" — 2" on either side of the phase center, 
again consistent with rotation. The largest velocity gradient is along a position angle 115°, 
i.e., close to the direction of the elongation of the continuum emission. In addition to this 
regular behavior, there also appears to be highly blueshifted emission (in the — 600kms~^ 
channel) northwest of the center. We interpret the line profiles and channel maps as showing 
rapid rotation about a position very near the phase center. Taken together, the line and 
continuum data indicate that the most likely morphology of the source is a disk of molecular 
gas and dust rotating about a central AGN located at the position of the hard X-ray peak. 
The absolute astrometry renders unlikely but does not exclude an alternative configuration 
of two galaxies orbiting each other, in which the AGN coincides with either the red or the 
blue-shifted CO(3-2) emission peak. 

4. Discussion 

4.1. J02399 is a very massive galaxy 

In what follows, we assume that the molecular gas in J02399 is indeed distributed 
in a rotating disk. We have constructed ring/disk models with flat rotation curves to flt 
the continuum and line proflle data, including the effects of beam smearing and nonzero 
velocity dispersions. The data are well fit (Figures 3 and 5) by a disk model with a (radially 
constant) rotation velocity of 420±20 km s~^ for an assumed inclination oii — 90°. The broad 
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wings of the integrated line profile indicate that large local random motions (turbulence) are 
present in the disk. The fit requires a ratio of local velocity dispersion to rotation velocity 
of ~ 0.23, which is ~ 2 — 3 times larger than in the Milky Way disk but comparable to 
that seen in local ULIRG mergers (Downes & Solomon 1998). A ring-like distribution with 
a maximum gas density at r ~ 1" on the sky {R = 3.2 h^lkpc in the source plane) and 
Ai?(FWHM)/i? ~ 1 — 1.5 provides a better match to the flux ratio of ~ 2 between the 
"horns" and the central dip in the global spectrum than a filled, centrally-peaked disk. 
Downes & Solomon (1998) have found that similar rings on somewhat smaller scales are 
characteristic of the molecular gas distributions in local ULIRGs. The outer disk radius 
required by our model is r ~ 2" — 2.5" (i? ~ 6 — 8/i^ykpc), in good agreement with our 
direct fit to the 335 /im continuum source size. The enclosed dynamical mass inferred from 
this disk model is 

Here 7 < 1 is a dimensionless scale factor which depends on the disk/spheroid structure, 
and which for a thin disk in a massive spheroid is close to unity (Lequeux 1983). Within the 
r ~ 2.5" outer radius of the molecular emission, the dynamical mass is ~ 3.2 x 10^^ h^} Mq. 
If a merger model is used instead of a rotating disk, the enclosed mass is approximately a 
factor of two larger. 

The integrated CO(3-2) fiux from J02399 is 3.1 ± 0.4 Jykms"^ (±0.3 Jykms"^ for a 
10% uncertainty in the fiux scale), in excellent agreement with the value of 3.0Jykms~^ 
found by Frayer et al. (1998). With this integrated fiux and the assumption of moderately 
dense {nn^ ^ lO^cm"^), thermalized gas at temperature > 30K (see e.g., Combes, Maoli, 
& Omont 1999), the molecular (plus helium) gas mass can be estimated following Solomon 
et al. (1997): 



Mq ■ Vjykms-V V2.45/ V9O.8IGHZ/ V23.5/io^Gpc/ V3.808. 

(2) 

Here a is the conversion factor from CO luminosity to gas mass, which we assume to be 
~ 0.8 — 1.6 Mq (K km pc^)~^, as estimated for local ULIRGs (Solomon et al. 1997; Downes 
& Solomon 1998). The inferred molecular gas mass in J02399 is thus 6.0±2.4 x 10^° h^^ Mq, 
implying a gas mass fraction of roughly 10-30% within r = 2.5" which is similar to those 
observed in local ULIRGs (Downes & Solomon 1998). 

J02399 is clearly a very massive system, and is in particular a very massive baryonic 
system. Observational studies of local ellipticals (Keeton 2001; Boriello, Salucci, & Danese 
2002; cf. Loewenstein & Mushotzky 2002) suggest that typically ~ 30% of the total mass 
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within the optical effective radius will be dark; the interior regions of massive local disk 
galaxies have similarly low dark matter fractions (Combes 2002). By analogy with these local 
systems, J02399 should have a baryonic mass of > 2 x 10^^ sin~^ i /iq 7 Mq within S/iQ^kpc; 
^gas ~ 6 X 10^'^ h^j Mq then implies that a conservative lower limit to the stellar mass 
within the same radius is ~ 1.4 x 10^^ Mq (for i = 90°). The Kennicutt (1983) stellar initial 
mass function (IMF) models of Cole et al. (2001)^ imply that m* - 7 x 10^° for the 

galaxy stellar mass function in the local Universe. Based on the stellar content of its inner 
8 h^^ kpc, then, J02399 is already at least a 2m* system, on its way to becoming a > 3m* 
system as it turns the remainder of its molecular gas into stars. 

Our new high-resolution observations of J02399 provide a fairly reliable estimate of the 
dynamical mass of a massive high- 2; galaxy. The more easily available observations of the 
ionized gas kinematics are more prone to systematic bias and influences other than gravity 
(Pettini et al. 2001; Heckman et al. 2000). The presence of such a massive object at a 
redshift as high as ~ 3 is qualitatively surprising, given that CDM models make the general 
prediction that galaxy masses build up gradually as a function of time and redshift, such 
that the comoving density of massive objects at 2; ~ 2 to 3 is only a small fraction of their 
local density (e.g., Kauffmann & Chariot 1998; Baugh et al. 2002). However, given the small 
number of objects whose kinematics have so far been studied with millimeter interferometry, 
the question arises whether the detection of a few such massive objects is significant, or 
whether they are mere items of curiosity. 

Besides J02399, the SCLS has yielded another source whose high redshift has been 
confirmed by spatially resolved CO spectroscopy. J14011 at zco = 2.565 has an observed 
850 //m flux density (12mJy), a lensing magnification (2.5: Ivison et al. 2001), a bolometric 
luminosity (1.3 x 10^^ Hq^Lq), and a molecular gas mass (3.4 x 10^^ Mq) which are 
similar to those of J02399 (Prayer et al. 1999; Ivison et al. 2001). The CO line width 
of J14011 (PWHM Av = 190 km s-^: Prayer et al. 1999; Downes & Solomon 2002) is 
much narrower than that of J02399 (PWHM Av = llOOkms"^), however. Although no 
clear rotation pattern has as yet been detected and no inclination correction can be made 
on the basis of the available CO line and millimeter continuum maps (Prayer et al. 1999; 
Downes & Solomon 2002), the inclination must be less than 20° for the dynamical mass to be 
> 7 X 10^° h^^ Mq (i.e., at least twice the gas mass). If the complex UV source structure is 
taken as evidence of a merger morphology, the dynamical mass would be > 1.4x10^^ h^^MQ. 
Downes & Solomon (2002) propose a disk model with a rotation velocity of 250kms~^, an 
inclination of 50°, and an outer radius of 1.1"; for these parameters and magnification by a 



^The Kennicutt (1983) IMF is comparable in this context to a Salpeter (1955) IMF from miower = 1 Mq to 
mupper = 100 Mq. For a 0.1 - 100 M© Salpeter IMF, Cole et al. (2001) find m*{z = 0) = 1.4 x 10^^ h^J^MQ. 
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factor M — 2.5, J14011's dynamical mass is 6 x 10^^ h^^ Mq. Their further suggestion that 
J14011 undergoes an additional factor of ten magnification by a previously unrecognized 
foreground galaxy (and that M^yn is therefore a factor of ten lower) rests on the assumption 
that the filling factor of CO emission within the l.l"-radius structure is unity. Of the seven 
background SCLS sources with intrinsic 5*850 > 4mJy (before lensing magnification), then, 
the two with known redshifts > 2 both have likely baryonic masses in the > 10^^ Mq bin. 
These source parameters are not atypical at high redshift: in a growing number of high-^ 
QSOs and radio galaxies, millimeter CO interferometry reveals comparably large molecular 
gas masses (in several cases already > lO^-*^ /i^y Mq) and a variety of hnewidths (Guilloteau 
et al. 1999; Papadopoulos et al. 2000; Cox et al. 2002, and references therein). 



4.2. Cosmic volume densities of massive galeixies as a function of redshift 

For our cosmology, the observed area covered by the SCLS (3.6 x 10~^sr) and the 
1 < z < 5 range sampled by the 850 //m observations correspond to a comoving volume 
of 5.3 X 10^/i^^Mpc^ The effective comoving volume density corresponding to one source 
like J02399 is thus $(M > 10^^ M©) = lO'^-^ hl^Mpc'^ . Including J14011, this becomes 
lO^^'^ /iq 7 Mpc~^; including the four remaining SCLS sources with intrinsic 5*850 ^ 4 mJy 
and unknown redshifts (on the assumption that they are also at 2; > 1 and comparably 
massive) boosts $ to lO^'^'^"'"^ /iq y Mpc^^. Although the latter step is the most uncertain, 
evidence from the weakness of radio and optical/UV counterparts for submillimeter sources 
suggests that these additional four sources are more likely to lie at high than at low redshifts 
(Dannerbauer et al. 2002; Ivison et al. 2002). Regardless of whether we assign one, two, or 
six SCLS sources to the > 10^^ Mq baryonic mass bin, however, we must adjust $ to reflect 
the fact that not all massive objects will be detectable by SCUBA as ~ lOmJy sources at all 
times. For a Miller-Scalo (1979) IMF, or a 1 - 100 M© Salpeter (1955) IMF, the bolometric 
(infrared) luminosity and star formation rate are related by (Kennicutt 1998) 

SFR Lit} , , 

r — (3) 

Moyr-i 1.2 x lO^o L© ^ ' 

such that the star formation rate in J02399 (for an assumed AGN contribution to Ljr of 
50%) is about 500 Mq yr"-*^. Turning a total baryonic mass of ~ 2 x 10^-*^ Mq into stars at this 
rate will require a total of ~ 4 x lO^yr. Since the time elapsed over the range 1 < 2; < 5 is 
~ 4.6 Gyr, correction for the finite starburst lifetime will increase the total $(M > 10^^ Mq) 
by at least a factor of ~ 11.5 over our estimates above. The correction factor can be higher 
still if any of the past star formation in J02399 occurred in a quiescent mode, if the current 
burst will be terminated prematurely by strong negative feedback to the interstellar medium 
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(e.g., Genzel & Cesarsky 2000), or if the dark matter fraction within R — 8 /iq.t kpc is greater 
than 30%. 

We now consider the evolution of $(M > 10^^ Mq) with redshift. In the upper panel of 
Figure 6, we have plotted at z ~ 2 — 3 the comoving volume densities derived for two and six 
SCLS sources having baryonic masses > 10^^ M©; upward arrows indicate a factor of 11.5 
starburst lifetime correction. At z ~ 1, we have plotted the volume densities of iC-selected 
systems with stellar masses > lO^-*^ h^^ Mq as inferred from a population synthesis analysis 
of Munich Near-lnfrarcd Cluster Survey (MUNICS) sources by Drory et al. (2002). Wc also 
show the volume densities of actively star-forming galaxies (S'15 > 300 /uJy) detected in deep 
ISOCAM surveys (Elbaz et al. 1999; Genzel & Cesarsky 2000). Rigopoulou et al. (2002) 
have obtained Ha rotation curves for three of these systems and find dynamical masses 
> 10^^ h^} Mq. Assuming that all ISOGAM galaxies with 15 //m flux densities > 300 //m in 
the redshift range 0.2 < 2; < 1.3 (sampled by the spectroscopy of Rigopoulou et al. 2002) 
have such masses, we derive a volume density of 2.4 x 10~^ /iq 7 Mpc~^. At z = 0, we show 
the comoving volume density of systems with stellar masses > 10^^ h^'^ Mq as obtained from 
modelling optical and near-IR 2dF and 2MASS photometry (Cole et al. 2001). 

For comparison with the observational data points, we also plot curves showing the pre- 
dictions of two different semi-analytic models for galaxy evolution within a AGDM structure 
formation cosmogony. For a representative "Munich" model, we count the number of galax- 
ies with total stellar+gas mass > lO^"*^ Mq in the catalogues of Kauffmann et al. (1999) and 
divide by the simulation volume of (200 /iq 7 Mpc)"^; the comoving densities are marginally 
higher than those we would derive from counting only systems with stellar mass > 10^^ Mq. 
For a representative "Durham" model, we use a curve for stellar mass > 10^^ Mq (C. Baugh, 
private communication) derived from a recent variation on the reference model of Cole et 
al. (2000). Although both Munich and Durham models adopt similar input parameters for 
their derivations of dark matter halo distributions {flu — 0.3, fl^ = 0.7, /10.7 = 1; cts — 0.9, 
and r ~ 0.2), they treat baryons differently Kauffmann et al. (1999) adopt a higher cosmic 
baryon density (f2b = 0.045 vs. 0.02), and the two groups use different recipes for star forma- 
tion, feedback, and other physical processes. As a result, the models' predictions for the rate 
of assembly of systems with baryonic masses > Mmin = 10^^ Mq differ substantially at high 
redshifts, as previously noted by Benson, Ellis, & Menanteau (2002)^°. Consideration of ob- 
servational data (note that Mmin scales differently with /iq.t for the different sets of sources) 
begins to make it possible to discriminate between the models. We find that the volume 
densities of massive SCUBA galaxies are clearly above the predictions of the Durham model. 



See in particular Figure 1 of astro-ph/01 10387. 
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and are significantly above the predictions of the Munich model if a reasonable correction for 
a finite starburst lifetime is applied. At ^ ~ 1, the data are still above the Durham model, 
but in good agreement with the Munich model. Both the Durham and the Munich models 
account well for the local volume densities of massive galaxies. 



J02399 is a type 2 QSO. It is thus of interest to compare our estimated volume densities 

for the brightest submillimctcr galaxies to those of high- 2; QSOs. Converting the UV/optical 
and X-ray luminosities of type 1 QSOs to bolometric luminosities with the Elvis et al. (1994) 
correction factors and assuming that the central black hole contains about 0.15% of the stellar 
spheroid mass (Ho 1999; Kormendy & Gebhardt 2001), the QSO host mass is related to the 
optical/X-ray luminosity by 



Here /j, is the bolometric correction for the band luminosity (/x ~ 30 and /opt ~ 14: 
Elvis et al. 1994), and r^Edd is the radiative efficiency relative to the Eddington rate.^^ For 
^Edd ~ 0.5, host masses > 10^^ /ig 7 M© correspond to X-ray luminosities > 3 x 10^^ /ig 7 erg s~^ 
and absolute B, V, and R magnitudes < —24.7. From the X-ray work of Miyaji, Hasinger, 
& Schmidt (2001) and the optical work of Hawkins & Veron (1996), we then estimate the 
density of type 1 AGN with Mhost > 10^^ h^^ Mq in the range 1.5 < ;2 < 3 to be ~ 
10-6±o.5;^3^Mpc~l The QSO lifetime tqso may be a factor of three shorter than tstarburst 
(Nusser & Silk 1993), and the ratio of type 2 to type 1 QSOs is p > 4 (GiUi, Salvati, & 
Hasinger 2001). Correcting for tstarburst/^QSO x p ~ 12, the comoving density of massive QSO 
hosts is ~ iQ-5±o.5 ^3 ^ Mpc~^, in remarkably- and probably fortuitously- good agreement 
with the values estimated for the bright SCUBA sources listed above. This agreement is 
consistent with the general assumption that the evolution of the most massive black holes 
and their host galaxies is connected through common formation at high redshift. 



Note that this relationship predicts a stellar spheroid mass for J02399 of > 2.4 x 10" Hq^Mq if 50% of 
its Lib, comes from a buried type 2 QSO, rjEdd ~ 0.5, and /ir > 1. 



4.3. Comparison to QSOs 




(4) 
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4.4. Is the discrepancy between observed and predicted volume densities 

mass- dependent ? 

CDM models do very well in accounting for the large-scale structure, luminosity func- 
tions, and color distributions of galaxies a.t z = (e.g., Kauffmann et al. 1999; Cole ct al. 
2000; papers in Guiderdoni et al. 2001). They predict the continuous buildup of galaxy 
masses from small to large through the basic process of hierarchical merging. If our ob- 
servations indicate that these models underpredict the high-z volume densities of massive 
galaxies, the next question is whether this shortfall originates in the predicted evolution 
of the underlying ACDM halo distributions, in some general feature of the semi-analytic 
models which build on them, or in some specific aspect of how the semi-analytic recipes 
treat baryonic processes in the most massive systems. Consideration of the halo catalogues 
of Kauffmann et al. (1999) offers no obvious evidence that the parameters of the ACDM 
simulations per se are at fault. If we multiply each halo mass by f2b (^m — ^b)~^ to estimate 
the associated baryonic mass, then even for the lower Qb — 0.02 preferred by the Durham 
group, the comoving volume densities of haloes with available baryonic mass > 10^^ M© 
comfortably exceed the volume densities of galaxies with observed baryonic mass > 10^^ Mq 
at every redshift. This statement remains true even if the point plotted for six SCLS sources 
in the upper panel of Figure 6 gains a full factor of ~ 11.5 correction for finite starburst 
lifetimes. 

To evaluate the performance of the semi-analytic models in a lower mass regime, we 
compare in the lower panel of Figure 6 the model predictions with the comoving volume 
densities of 2; ~ 3 Lyman break galaxies (LBGs) . These systems have average stellar masses 
near but somewhat below the knee of the Schechter function (M ~ 10^^ h^j Mq ~ 0.14m*), 
according to the recent stellar synthesis analyses of Papovich, Dickinson, & Ferguson (2001) 
and Shapley et al. (2001). Again, we compare these measurements to the Durham and 
Munich model predictions. In this lower mass bin, the semi-analytic models still fall below 
(but now closer to) the observed volume densities of LBGs. In addition, we have plotted 
the model predictions of Somerville, Primack, & Faber (2001), who include the effects of 
extinction on the predicted volume densities for R < 25.5 LBGs and consider two extreme 
star formation histories- one with a constant star formation timescale, and one including star 
formation bursts triggered by galaxy collisions and mergers. Clearly the low star formation 
rates predicted by the constant model fail by a large factor in accounting for the data, while 
the burst models fit the data very well. The Munich and Durham models (which include 
some acceleration of the star formation rate at high redshift) lie between these two extremes. 
It appears that the recipes for star formation and feedback used in the semi-analytic models 
do reasonably well in accounting for sub-m* galaxies, once the influence of merger-induced 
bursts is included. 
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Our discussion thus suggests that the discrepancy between models and data may increase 
with galaxy mass. Systems with masses a few times m* have comoving number densities 
at z ~ 1 which are comparable to those in the local Universe {z = 0). At z = 2 — 3, the 
same kinds of objects still have comoving volume densities which are > 10% of their z = 
values. These findings appear to contradict current hierarchical CDM models, which predict 
a large decrease of the abundance of massive galaxies at z > 1 due to an extended phase 
of merger activity (e.g., Kauffmann & Chariot 1998). The data presented here indicate 
that this picture may have to be revised with respect to the rate at which massive systems 
(and bulges) are assembled at fairly high redshift, possibly through adjustment of and 
its coupling to very efficient star formation in major mergers. More work is needed before 
a firmer statement can be made about this crucial, "baryonic mass assembly" test of the 
current standard cosmogony. Both observational and theoretical estimates are still very 
uncertain, with the former depending on small samples and large lifetime corrections, and 
the latter on ad hoc input recipes for star formation and feedback. 

We thank IRAM staff members, particularly Roberto Neri, for their help in coordinating 
and conducting the observations for this project, and Rob Ivison for sharing information 
about the SCLS lensing model. We are especially grateful to Len Cowie for providing us 
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thank Carlton Baugh for providing the most current Durham model predictions, Mark Bautz 
for clarifying the Chandra astrometry, and Chung-Pei Ma, Chris McKee, and Ian Small 
for very useful comments. This research has made use of the NASA/IPAC Extragalactic 
Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute 
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Table 1. Positions and uncertainties for J02399 in different wavebands. 



Band 




R.A. 


Dec. 


la astrometric 




(rest-frame) 


Feature 


(J2000) 


(J2000) 


uncertainty 


Reference 


335 fim 


centroid 


02:39:51.87 


-01:35:58.8 


±0.6" 


1 


1685 A 


peak (LI) 


02:39:51.845 


-01:35:58.2 


±0.3" 


2" 


0.4-1.8 keV 


peak 


02:39:51.81 


-01:35:58.1 


±0.6" 


3^ 



^Based on an astrometric solution of USN0-A2 stars relative to VLA sources in the 
Abell 370 field. 

'^Includes a (AR.A., ADec) = (—1.6", ±0.5") overall "boresight" correction between 
Chandra and APM optical reference frames. 



References. — (1) this paper (2) L. Cowie, private communication (3) Bautz et al. 
(2000); M. Bautz, private communication 
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Fig. 1. — Contours of 1.27 mm (rest- frame 335 /im) continuum emission from J02399 mapped 
with the PdBI, superposed on a greyscale image of the i?-band (rest-frame UV) contin- 
uum obtained at Keck (L. Cowie, private communication). Contours are multiples of 
0.4mJybeam~^, for a synthesized beam of 1.8" x 1.4" at position angle 42° east of north 
(shown at lower left). Arrows indicate rest-UV peaks LI and L2; cross denotes the position 
and la positional uncertainty (including both X-ray and millimeter astrometric uncertain- 
ties) of the hard X-ray source (Bautz et al. 2000; M. Bautz, private communication). The 
absolute uncertainty of the PdBI astrometry is ±0.6"; the i?-band astrometry is based on 
a comparison of USN0-A2 stars to VLA sources in the Abell 370 field and has a relative 
uncertainty of ±0.3"; the X-ray astrometry has an uncertainty of ±0.6". 

Fig. 2. — Overlay of continuum (contours, as in Figure 1) and CO (3-2) velocity-integrated 
line emission (greyscale). The CO(3-2) map has a synthesized beam of 5.2" x 2.4" at position 
angle 33° east of north. 

Fig. 3. — CO(3-2) spectrum, smoothed to 60kms^^ resolution and integrated over the 
central ~ 5" of the source. Per-channcl rms is indicated at lower left. Velocities arc with 
respect to the assumed redshift zco = 2.808 (Frayer et al. 1998). The dashed curve is a model 
spectrum, computed for a rotating disk with constant rotation velocity v^ot = 420kms~^, 
local FWHM random motions of 220kms~^, inclination i — 90°, and a Gaussian radial 
intensity distribution which peaks at 1.0" and has radial FWHM 1.5" (see §4.1). The FWHM 
vertical thickness of the disk/ring is assumed to be 0.5". The model takes into account spatial 
and velocity smearing. 

Fig. 4. — CO(3-2) channel maps in steps of 200kms~^ with 0.35mJybeam~^ contour spac- 
ing. Velocities are with respect to the assumed redshift zco — 2.808 (Frayer et al. 1998). 
The synthesized beam (shown at lower left) is 5.2" x 2.4" at position angle 33° east of north. 

Fig. 5. — Position- velocity diagram of observed CO(3-2) emission, at 60kms^^ resolution 
and position angle 115° (white contours are multiples of 0.55 mJy beam"^), overlaid on the 
predictions of the rotating disk model (greyscale) described in §4.1 and in the caption for 
Figure 3. 

Fig. 6. — Cosmic densities of galaxies of (baryonic/stellar) mass > Mmin for the fiat Q\ = 0.7, 
Hq = 70 hoj kms~^ Mpc~^ cosmology used throughout. Upper panel: M^i^ = 1 x 10^^ Mq. 
The filled dark circles (and la error bars) are estimated from the SCLS galaxies J02399 and 
J14011 (lower), and from the six SCLS galaxies with intrinsic Ss5o > 4 mJy which may lie at 
z > 1 (upper). Mass estimates are from the CO kinematics in J02399 and J14011 (this paper; 
Ivison et al. 2001; Downes & Solomon 2002). The dashed upward arrows indicate correction 
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factors for finite starburst lifetimes, as discussed in §4.2. The light filled square denotes the 
density of 5*15 > 300 yuJy ISOCAM galaxies in the range z = 0.2 — 1.3, for three of which 
Rigopoulou et al. (2002) have determined large dynamical masses from Ha rotation curves. 
The open rectangles with crosses are source densities for massive early-type galaxies from the 
MUNICS K-hand work of Drory et al. (2002), where stellar masses are estimated from stellar 
synthesis models. The filled triangle and bar denote the z = results obtained from stellar 
synthesis modelling of 2dF/2MASS data by Cole et al. (2000). For comparison, we show 
two theoretical predictions from semi-analytic models in CDM cosmogonies: the continuous 
curve is from a "Durham" model (C. Baugh, private communication; see also Cole et al. 
2001; Baugh et al. 2002), while the dashed curve is from a "Munich" model (Kauffmann et 
al. 1999). In all cases, the stellar masses refer to a Miller-Scalo (1979) IMF, or a Salpeter 
(1955) IMF between 1 and 100 M©. Lower panel: M„ii„ = 10^° Mq. Estimates of the 
density of 2; ~ 3 Lyman break galaxies with M > Mmin (as determined from stellar synthesis 
modelling) arc denoted by a filled rectangle (Shapley et al. 2001) and downward-pointing 
light triangle (Papovich ct al. 2001). All other notation is as in the upper panel. Thin and 
thick short-dashed curves are models from the work of Somerville et al. (2001), which take 
into account dust extinction and represent constant and merger- induced ( "coUisional" ) star 
formation histories, respectively. 
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